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The  reflection  and  scattering  properties  of  an  inhomogeneous  poroelastic  medium  were  studied 
via  numerical  simulation.  The  inhomogeneous  medium  was  modeled  as  an  ensemble  average  of 
randomly  layered  poroelastic  material.  Each  layer  represented  a  granular  material  of  a  particular 
grain  size.  The  thickness  of  each  layer  was  related  to  the  associated  grain  size  and  porosity  by  a 
conservation  of  mass  relationship.  Lateral  variations  in  grain  size  were  approximated  by 
performing  a  coherent  ensemble  average  of  results  from  several  realizations  of  the  randomly 
stratified  medium.  Poroelastic  medium  parameters  were  chosen  to  represent  water-saturated 
sand.  The  mean  and  standard  deviation  of  the  grain  size  distribution  were  chosen  to  match 
existing  experimental  data  so  that  the  model  could  be  tested.  Specifically,  the  inhomogeneous 
medium  was  modeled  as  bounded  by  a  homogeneous  water  half-space  on  the  source  side,  and  a 
homogeneous  poroelastic  half-space  of  equivalent  average  porosity  on  the  other  side.  Reflected 
signals  were  computed  for  500  kHz  and  1  MHz  normally  incident  plane  waves.  Coherent  and 
random  components  of  the  reflected  signal  were  calculated.  The  coherent  part  was  directly 
related  to  the  reflection  coefficient.  The  random  component  was  related  to  the  scattering  strength 
of  the  medium.  It  was  found  to  increase  with  the  mean  grain  size  diameter,  consistent  with 
previous  experimental  results. 
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PREFACE 


This  is  the  final  report  on  work  that  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin  (ARL:UT),  was  tasked  to  perform  under  Grant 
N0001 4-94-1 -0438,  entitled  "Multiple  Scatter  Theory  of  Ocean  Sediments." 
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1.  OBJECTIVE 


The  goal  of  this  study  is  to  investigate  the  extent  to  which  sediment 
granularity  can  account  for  the  attenuation  and  scattering  that  has  been 
observed  in  ocean  sediments,  particularly  water-saturated  sand.  Typically,  a 
couple  of  devices  are  employed  to  match  theory  to  experimental  measurements 
of  attenuation  and  scattering;  the  attenuation  is  accounted  for  in  terms  of  a 
complex  bulk  modulus,  and  the  scattering  is  modeled  as  a  random  field  of  point 
scatterers,  in  which  the  scattering  strength  is  adjusted  to  match  the  measured 
data.  Both  devices  have  very  little  physical  basis.  Our  hypothesis  is  that  both 
attenuation  and  scattering  might  be  explained  in  terms  of  acoustic  interaction 
with  the  sediment  grains.  The  viscoelastic  theory  of  acoustic  propagation  in  a 
solid  is  not  a  suitable  starting  point  since  it  does  not  possess  a  mechanism  that 
can  account  for  the  interaction  between  pore  fluid  and  solid  particles.  Our 
starting  point  is  Biot's  theory  of  acoustic  propagation  in  a  poroelastic  medium.  It 
contains  the  basic  mechanisms  of  acoustic  interaction  between  the  solid  matrix 
and  pore  liquid  as  far  as  forward  propagation  is  concerned,  but  it  does  not  have 
any  mechanism  to  account  for  scattering  and  any  associated  losses  due  to 
granularity.  Our  intention  is  to  extend  Biot's  theory  to  inciude  the  effects  of 
granularity.  In  this  initial  study,  we  will  use  a  brute  force,  numerical  simulation 
approach  to  obtain  results  quickly.  The  results  of  this  initial  study  were 
compared  with  experimental  data  to  test  the  feasibility  of  our  hypothesis.  The 
comparison,  as  it  turned  out,  was  very  encouraging,  and  provides  a  solid 
foundation  for  proceeding  to  the  next  phase  of  the  investigation,  which  will  be 
directed  towards  obtaining  analytical  relationships. 
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2.  INTRODUCTION 


Current  models  of  acoustic  bottom  backscatter  from  sandy  sediments  are 
based  on  composite  roughness  and  volume  scattering  effects  from  sediments 
modeled  as  fluids.^'®  However,  recent  bottom  penetration  experiments  by 
Chotiros"^'®  have  shown  that  Biot's  theory,^’®  based  on  the  propagation  of 
sound  through  a  fluid-filled  porous  solid  matrix,  better  explains  the  conduction 
of  acoustic  energy  into  the  sediment.  This  conduction  process  is  an  important 
factor  in  any  model  of  scattering  by  sediment  volume,  making  it  plausible  to 
apply  Biot's  theory  to  sediment  scattering.  Biot's  theory  predicts  acoustic 
propagation  through  a  poroelastic  medium  as  a  coupled  wave  motion  within  the 
solid  and  the  pore  fluid.  The  medium  is  modeled  as  a  solid  with  tubular  pores, 
as  illustrated  in  Fig.  2.1.  The  propagating  wave,  in  the  direction  of  the  pore 
tubes,  can  be  decomposed  into  three  components  -  a  fast  and  a  slow 
compression  wave,  and  a  shear  wave.  The  equations  of  motion,  in  a  form  used 
by  Stern,  Bedford,  and  Millwater®  (correcting  for  typographical  errors),  are 

+  (//-^)V(V-u)  -  C  V(V- w)  =  pii  -  p^w  ,  (2.1) 

CV(V-u)  -  MV(V-w)  =  p^ii  - -^w  - ,  (2.2) 

p  K 

where  u  is  the  displacement  vector  of  the  solid  frame,  w  is  the  negative  porosity 
times  the  displacement  vector  of  the  pore  fluid  relative  to  the  solid  frame,  p  is 
the  porosity  of  the  solid  frame,  pf  is  the  mass  density  of  the  pore  fluid,  p  is  the 
mass  density  of  the  saturated  sediment,  p  is  the  shear  modulus  of  the  solid 
frame,  c  is  the  virtual  mass  coefficient  of  fluid  motion,  t]  is  the  viscosity  of  the 
fluid,  and  k  is  the  permeability  of  the  solid  frame.  C,  H,  and  M  are  constitutive 
coefficients  depending  on  p,  p,  and  the  bulk  moduli  of  the  pore  fluid,  grain 
material,  and  the  saturated  sediment.  F  is  a  frequency-dependent  dynamic 
correction  term  describing  the  frictionai  force  due  to  the  relative  motion  of  the 
solid  and  fluid,  as  given  by  Biot.^’®  The  relationships  between  the  parameters 
of  Eqs.  (2.1)  and  (2.2)  are  given  id  Appendix  A.  In  order  to  match  the  calculated 
acoustic  wave  attenuations  to  measured  values,  imaginary  terms  must  be 
added  to  the  frame  bulk  and  shear  moduli.  Furthermore,  the  theory  has  no 
mechanism  for  scattering  because  the  pores  are  modeled  as  parallel  smooth 
walled  tubes  that  are  perfectly  aligned  with  the  direction  of  wave  propagation. 
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When  u  and  w  are  written  in  terms  of  vector  and  scalar  potentials,  the 
above  two  equations  separate  into  four  coupled  vector  partial  differential 
equations.  For  homogeneous  media,  the  vector  and  scalar  potentials  have 
plane  wave  solutions.  At  a  fluid/poroelastic  interface,  there  are  four  boundary 
conditions  imposed  on  these  solutions  --  continuity  of  fluid  pressure,  continuity 
of  shear  traction,  continuity  of  normal  traction,  and  continuity  of  normal  fluid 
displacement.  At  a  poroelastic/poroelastic  interface  there  are  two  additional 
boundary  conditions  --  continuity  of  tangential  solid  displacement  and  continuity 
of  normal  solid  displacement.  The  mathematical  formulation  of  the  differential 
equations  and  boundary  conditions  governing  the  vector  and  scalar  potentials 
are  clearly  illustrated  by  Stern,  Bedford,  and  Millwater,®  and  are  outlined  in 
Appendix  A. 
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3.  MODEL 


We  now  present  a  multiple  scatter  model  of  the  ocean  sediment.  In  this 
model,  the  ocean  sediment  is  represented  as  a  finite  granular  medium  bounded 
above  by  a  homogeneous  fluid  half-space  (water)  and  below  by  a 
homogeneous  poroeiastic  half-space.  The  poroelastic  half-space  has 
properties  given  in  Table  3.1,  which  were  recently  determined  by  Chotiros^  for 
sound  propagating  through  water-saturated  sand.  Likewise,  the  fluid  half-space 
has  properties  denoted  by  the  fluid  properties  in  Table  3.1.  Before  continuing, 
the  terms  “granular"  and  "poroelastic"  need  some  clarification.  In  the  context  of 
this  study,  the  term  "poroelastic"  refers  to  a  medium  that  follows  Biot's  theory  of 
acoustic  propagation  as  described  by  Eqs.  (3.1)  and  (3.2),  in  which  the  pores 
are  modeled  as  parallel  tubes  as  described  in  Fig.  2.1.  The  term  "granular" 
refers  to  a  poroelastic  medium  that  is  consistent  with  Biot's  theory,  and  which,  in 
addition,  the  pore  diameter  may  vary,  thus  giving  rise  to  scattering  mechanisms 
associated  with  the  granularity. 

A  randomly  packed  granular  medium,  such  as  the  ocean  sediment,  is 
essentially  a  three-dimensional  medium  that  will  reflect,  transmit  and,  most 
importantly,  scatter  an  incident  acoustic  wave.  Ideally,  a  fully  three-dimensional 
mathematical  model  is  required,  but  it  would  be  very  expensive  to  develop  and 
compute.  However,  it  is  noted  that  the  inhomogeneities  in  ocean  sediments, 
particularly  sand  beds,  are  structured  as  horizontal  lenses,  as  illustrated  in 
Fig.  3.1(a).  We  have  chosen  to  model  the  lenses  as  an  ensemble  of  random, 
horizontally  layered,  poroelastic  media,  as  illustrated  in  Fig.  3.1(b),  whose 
acoustical  response,  when  summed,  will  approximate  that  of  a  fully  three- 
dimensional  model.  This  model  ignores  the  edge  effects  at  the  periphery  of  the 
lenses,  which  are  expected  to  be  small.  The  approach  is  attractive  because  a 
horizontally  layered  medium,  being  essentially  a  one-dimensional  problem,  is 
quite  simple  to  compute,  as  shown  in  Appendix  B.  Even  an  ensemble  of 
several  random  realizations  is  still  quite  economical  computationally. 

Each  realization  of  a  randomly  layered  poroelastic  medium  consists  of  a 
series  of  poroelastic  layers  of  differing  properties,  sandwiched  between  an 
upper  fluid  half-space  and  a  lower  uniform  poroelastic  half-space,  as  illustrated 
in  Fig.  3.2.  Within  the  randomly  layered  medium,  the  layer  corresponds  to 
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Table  3.1 

Material  properties  of  the  homogeneous  Biot  half-space. 


Fluid  viscosity  (77) 

Fluid  mass  density  {pf) 

Fluid  bulk  modulus  (Kf) 

I.OOx  lO'^kg/m-s 

1000  kg/m3 

2.25x109  Pa 

Grain  mass  density  (pg) 

Grain  bulk  modulus  (/C) 

2650  kg/m3 

7.00  X  109  Pa 

Frame  shear  modulus  {po) 
Frame  bulk  modulus  {Kbo) 

Frame  porosity  (jS) 

Frame  shear  log  decrement  (5s) 
Frame  bulk  log  decrement  (5) 

2.61  X  107  Pa 

5.30x109  Pa 

0.36 

0.15 

0.15 
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Random  scattered  field 


Incident  plane  wave 


(a)  Random  scattering  of  a  plane  wave  from  a  3-D  granular  structure. 


Plane  wave  responses 
from  random  realizations 


(b)  Summation  of  scattered  signals  from  ensemble  of  1-D  structures. 


Figure  3.1 

Simulation  of  plane  wave  scattering  from  a  3-D 
granular  structure  with  an  ensemble  of  1-D  structures. 


ARL - UT 
AS-95-1101 
DJY/RFO 
12-95/svr 


9 


Upper  fluid  half-space 


■rVr~M 


Tubular 
pores  of 

variable  diameter 
in  a  solid 
matrix 


Lower  half-space 
of  tubular  pores  of 
constant  diameter 
in  a  solid  matrix 


Figure  3.2 

A  randomly  layered  Biot  medium,  bounded  above  by  a 
homogeneous  fluid  half-space  and  below  by  a 
homogeneous  Biot  half-space. 
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a  mono-layer  of  granular  material  of  uniform,  but  randomly  chosen,  grain 
diameter,  a;.  For  a  given  grain  size,  greater  layer  thickness  implies  greater  fluid 
volume  within  the  layer,  and  hence  greater  porosity  for  the  layer.  Therefore,  a 
certain  correlation  between  the  thickness  hyand  porosity  j3yof  the  itii  iayer  within 
the  granular  medium  is  postulated.  The  layer  porosity  is  simply  the  ratio  of  the 
fluid  volume  to  the  total  volume  of  the  layer. 


The  average  layer  thickness  is  chosen  so  that  average  porosity  of  the 
granular  medium  is  the  same  as  that  of  the  homogeneous  poroeiastic  half¬ 
space,  as  given  in  Table  3.1. 

The  properties  that  are  dependent  on  porosity  and  grain  diameter  are 
computed  according  to  the  method  described  by  Chotiros.^  Other  material 
properties  of  each  layer  within  the  granular  medium  match  those  of  the 
homogeneous  Biot  half-space.  Further  details  are  given  in  Appendix  A.  Of 
particular  significance,  within  the  randomiy  layered  poroeiastic  medium,  the  log 
decrements  are  set  to  zero  (i.e.,  the  imaginary  parts  of  the  bulk  and  shear 
moduli  are  set  to  zero),  in  the  expectation  that  the  acoustic  interactions  with  the 
pore  size  variations  wili  provide  the  necessary  loss  and  scattering  mechanism 
that  will  properly  account  for  the  experimentally  observed  attenuation  and 
scattering  phenomena. 
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4.  SIMULATION  PROCEDURE 


The  purpose  of  this  simulation  was  to  verify  the  feasibility  of  the  model 
described  above  by  comparing  the  model  predictions  of  reflection  loss  and 
scattering  strength  with  experimental  data.  The  simulation  parameters  were 
chosen  to  allow  comparison  of  the  results  with  the  experimental  work  of  Nolle^° 
and  Mifsud,^^  who  made  measurements  at  500  kHz  and  1  MHz.  The  mean  and 
standard  deviations  of  the  grain  sizes  of  the  sand  samples  used  in  their 
experiments  are  shown  in  Table  4.1.  For  simplicity,  only  the  case  of  normal 
incidence  was  considered. 

It  was  assumed  that  the  grain  size  distribution  is  a  log-normal  function. 
Hence,  the  distribution  is  completely  specified  by  two  free  parameters,  its  mean 
and  standard  deviations.  For  simulation  purposes,  a  number  of  mean  grain 
sizes  between  120  pm  and  640  pm  were  chosen.  The  dimensionless  standard 
deviation  was  then  determined  by  linear  interpolation  of  the  corresponding 
experimental  values  in  Table  4.1. 


Table  4.1 

The  mean  and  dimensionless  standard  deviations  for  the  grain  size 
distributions  of  the  granuiar  medium.  (The  dimensioniess  standard 
deviations  for  intermediate  grain  diameters  are  determined  by 

iinear  interpoiation.) 


Mean  diameter 

Dimensionless  standard 

(urn) 

deviation 

120 

0.28 

170 

0.17 

400 

0.15 

640 

0.11 

In  the  absence  of  any  further  information,  the  dimensionless  layer 
thickness  was  arbitrarily  chosen  to  have  the  same  statistical  distribution  as  the 
grain  size,  but  the  grain  size  and  layer  thickness  were  treated  as  independent 
variables.  The  mean  of  the  thickness  distribution  was  chosen  so  that  the  mean 
porosity  of  a  many-layered  sediment,  as  determined  from  Eq.  (4.1),  was  the 
same  as  that  of  the  homogeneous  poroelastic  half-space,  i.e.,  0.36.  Thus, 
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IIRII 


Total  thickness  -  mean  grain  diameter 


Figure  4.1 

An  example  showing  the  variation  of  the  reflection  amplitude  with 
the  total  sediment  thickness  (mean  grain  diameter  =  120  )xm). 
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1  -  0.36 


0.64 


(4.1) 


Given  the  grain  size  and  layer  thickness  distributions,  the  layered  Biot 
medium  was  generated  layer  by  layer.  The  grain  size  and  thickness  of  each 
layer  were  independently  chosen  at  random,  consistent  with  the  given 
distribution  functions.  The  reflection  loss  was  computed  as  a  function  of  the 
number  of  layers,  using  the  procedure  given  in  Appendix  A.  One  realization  of 
such  a  plot  is  shown  in  Fig.  4.1.  When  the  number  of  layers  was  small,  the 
reflection  amplitude  varied  wildly  with  the  addition  of  each  new  layer.  As  the 
number  of  layers  increased,  the  amplitude  of  the  reflected  signal  appeared  to 
converge  about  a  mean  value. 

To  account  for  lateral  variations  in  a  realistic  sandy  sediment,  the 
reflected  signal  was  computed  from  an  ensemble  of  random  realizations  of  the 
same  grain  size  and  layer  thickness  distributions.  The  ensemble  average  of  the 
reflected  signal  was  used  to  compute  the  reflection  coefficient.  It  was  implicitly 
assumed  that  there  is  no  correlation  between  the  lateral  depth  profile  variations. 
Typically,  25  random  realizations  were  generated  for  each  grain  size 
distribution.  The  reflection  coefficient  was  calculated  from  the  mean  value  of  the 
complex  reflected  signal,  i.e.,  its  coherent  component,  as  given  by  ll<R>ll^.  The 
scattering  strength  was  calculated  from  the  variation  in  reflected  amplitudes,  i.e., 
the  random  component,  as  given  by 

S  =  101og((||fi-(/()f))  =  101og((||sr)-pf)  .  (4,2) 

Due  to  the  one-dimensional  nature  of  these  simulations,  the  scattered 
energy  was  all  directed  back  to  the  source  of  the  incident  wave.  But, 
experimentally,  the  measured  scattering  at  a  water-sediment  interface  is  a 
three-dimensional  process  in  which  energy  is  scattered  within  a  hemispherical 
sector  of  2;c  steradians.  To  compensate  for  this  inherent  discrepancy,  the 
simulated  scattering  strength  was  offset  by  a  precalculated  amount  that 
depends  on  the  incident  frequency.  The  details  of  this  calculation  are  given  in 
Appendix  C. 
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5.  RESULTS 


5.1  REFLECTION  LOSS 

As  illustrated  in  Fig.  4.1,  the  norm  of  the  reflection  amplitude  varied  about 
a  mean  value  as  the  number  of  layers  in  the  sediment  increased.  The  values 
were  averaged  over  575  realizations  of  the  granular  sediment.  The  result  was 
estimated  to  be  0.268,  equivalent  to  -11.4  dB,  which  is  in  good  agreement  with 
the  measured  value  by  Nolle^®  and  Mifsud^^  of  -11  dB.  No  significant 
deviations  from  this  value  were  observed  between  the  frequencies  or  grain 
sizes  used  in  this  study. 


5.2  SCATTERING  STRENGTH 

For  incident  wave  frequencies  of  both  500  kHz  and  1  MHz,  scattering 
strengths  were  calculated  for  mean  grain  diameters  ranging  from  120  pm  to 
640  pm.  The  results  are  shown  superimposed  on  the  experimental  data  of 
Nolle  in  Figs.  5.1  and  5.2. 

At  500  kHz,  the  calculated  scattering  strength  at  normal  incidence  as  a 
function  of  grain  size  was  in  good  agreement  with  the  experimental  data  for 
mean  grain  diameters  up  to  about  530  pm.  Beyond  this,  the  calculated 
scattering  strength  dropped  sharply  below  the  experimental  values.  At  the  point 
of  divergence,  the  wavelength  of  the  Biot  fast  wave  was  6.4  mean  grain 
diameters,  or  about  four  mean  layer  thicknesses.  This  indicates  that  our  model 
is  only  valid  for  grain  sizes  significantly  less  than  the  acoustic  wavelength  within 
the  sediment. 

At  1  MHz,  the  calculated  scattering  strength  was  only  in  agreement  with 
the  experimental  data  for  mean  grain  diameters  up  to  about  210  pm.  At  this 
point  of  divergence,  the  wavelength  of  the  Biot  fast  wave  was  8.0  mean  grain 
diameters,  or  about  five  mean  layer  thicknesses.  Although  these  results  are 
consistent  with  the  previous  interpretation  that  the  model  is  valid  for  grain  sizes 
significantly  less  than  the  acoustic  wavelength,  the  range  of  experimental  data 
points  below  210  pm  was  too  small  to  establish  a  legitimate  trend. 
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Figure  5.1 

Calculated  backscattering  strength  as  a  function  of  mean  grain 
diameter  at  500  kHz  compared  with  Nolle's  experimental  data. 
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Figure  5.2 

Calculated  backscattering  strength  as  a  function  of  mean  grain 
diameter  at  1  MHz  compared  with  Nolle's  experimental  data. 


ARL - UT 
AS-95-1105 
DJY/RFO 
12-95/svr 


19 


This  page  intentionally  left  blank. 


20 


6.  CONCLUSIONS 


The  effects  of  granularity  have  been  added  to  a  Biot  model  of  ocean 
sediments  in  a  numerical  simulation  by  introducing  homogeneous  layers  of 
random  thicknesses  on  the  order  of  a  grain  diameter,  in  which  the  mean 
porosity  is  preserved.  Within  each  layer,  conservation  of  grain  volume  dictates 
a  linear  relationship  between  the  layer  thickness,  porosity,  and  grain  diameter. 

Lateral  variations  in  sediment  structure  were  simulated  by  coherently 
averaging  the  results  for  several  random  realizations  of  a  layered  poroelastic 
medium  with  given  grain  size  and  layer  thickness  distributions.  The  reflection 
loss  predicted  by  this  model  was  computed  from  the  coherent  component  of  the 
ensemble  average  of  the  reflected  signal,  and  it  was  found  to  be  in  good 
agreement  with  measured  values  by  Nolle.  The  scattering  strength  was 
computed  from  the  random  component  of  the  reflected  signal,  and  its  values 
were  found  to  be  in  agreement  with  experimental  data  for  mean  grain  sizes 
significantly  less  than  the  acoustic  wavelength  within  the  sediment.  However, 
the  simulated  scattering  strength  dropped  sharply  below  the  observed  values 
for  grain  sizes  near  or  greater  than  the  acoustic  wavelength.  This  indicates  that 
the  model  is  only  valid  for  incident  waves  that  yield  acoustic  wavelengths  near 
or  greater  than  ten  times  the  mean  grain  diameter  of  the  sediment. 

The  approach  that  we  have  taken  is  based  on  sound  physical  principles. 
It  is  a  significant  extension  of  Biot's  theory  of  acoustic  propagation  in  porous 
media,  and  gives  an  insight  to  the  processes  that  give  rise  to  reflection  and 
scattering  from  a  granular  medium  such  as  water  saturated  sand,  that  is  directly 
applicable  to  ocean  sediments.  The  above  agreement  between  model 
predictions  and  experimental  measurement  indicates  that  the  approach  is 
feasible.  Follow-on  work  will  include  computation  of  wave  attenuation  due  to 
granularity,  which  would  involve  calculation  of  transmission  coefficients  through 
the  granular  material.  Then,  we  will  attempt  to  derive  an  analytical  model,  based 
on  our  understanding  of  the  numerical  simulation  results. 
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APPENDIX  A 

PROBLEM  SETUP  AND  SOLUTION 
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Most  of  the  equations  used  are  the  same  as  those  given  by  Stern, 
Bedford,  and  Millwater  ®  There  are  a  few  typographical  errors  in  the  equations 
from  this  reference,  which  have  been  corrected.  The  corrected  equations  in  this 
appendix  are  Eqs.  (A.5),  (A.6),  (A.29),  (A. 45),  and  (A.46). 

The  ocean  sediment  is  represented  as  a  finite  granular  medium  bounded 
above  by  a  homogeneous  fluid  half-space  (water)  and  below  by  a 
homogeneous  Blot  half-space.  The  granular  medium  is  composed  of 
homogeneous  Biot  layers  of  differing  properties  and  thicknesses,  as  depicted  in 
Fig.  3.2. 

The  layers  of  the  finite  granular  medium  are  numbered  1  to  N,  from  top  to 
bottom.  The  homogeneous  Biot  half-space  Is  layer  number  N  +  ^,  and  the  fluid 
medium  is  number  0.  All  /V  +  1  Biot  layers  have  identical  values  for  the  fluid 
viscosity  (77),  fluid  mass  density  (pf),  fluid  bulk  modulus  (Kf),  grain  mass  density 
(pg),  grain  bulk  modulus  (Kr),  frame  shear  modulus  (po),  and  frame  bulk 
modulus  {Kbo).  The  fluid  medium  has  the  same  fluid  mass  density  {pf)  and  fluid 
bulk  modulus  {Kfj  as  the  Biot  media.  These  values  are  given  in  Table  3.1 . 

The  '0-  layer  of  the  granular  medium  corresponds  to  a  mono-layer  of 
granular  material  of  uniform,  but  randomly  chosen,  grain  size,  a/.  The  choice  of 
a  grain  size  distribution  is  the  sole  free  parameter  for  generating  a  granular 
sediment  of  given  total  thickness. 

In  this  model  there  is  a  correlation  between  the  thickness  hi  and  porosity 
Pi  of  the  '0  layer  within  the  granular  medium.  Each  layer  represents  a  mono- 
layer  of  granular  material  of  uniform  grain  size.  So,  for  a  given  grain  size, 
greater  layer  thickness  implies  greater  fluid  volume  within  the  layer,  and  hence 
greater  porosity  for  the  layer.  For  a  given  grain  size  distribution,  the 
dimensionless  layer  thickness  distribution  is  identical  to  that  of  the  grain  size 
distribution.  The  mean  of  the  layer  thickness  distribution  is  chosen  so  that 
average  porosity  of  a  many-layered  granular  medium  would  be  the  same  as  for 
the  homogeneous  Biot  half-space,  pn+i.  The  porosity  of  the  itli  layer  is  simply 
the  ratio  of  the  fluid  volume  to  the  total  volume  of  the  layer.  Due  to  the 
effectively  one-dimensional  character  of  the  simulations  in  this  present  study. 
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we  choose 


Pi  =  1 


i  =  l,2y...,N  , 


(A.1) 


Pn  +  i  —  0.36. 


(A.2) 


The  skeletal  frame's  bulk  and  shear  logarithmic  decrements  are  set 
equal  to  zero  for  the  first  N  layers,  and  set  equal  to  0.15  for  the  final  semi¬ 
infinite  layer. 


5i  =  8,i  =  0,  i  =  , 


(A.3) 


5n  +  \  —  SsN  +  i  0*15  , 


(A.4) 


Kbi  “■  1  / 


(A.5) 


(A.6) 


The  mass  density  of  the  saturated  sediment  is  simply 


Pi  =  (1  -  ft  dps  +  Pipf 


(A.7) 


The  constitutive  coefficients  given  in  Eqs.  (2.1)  and  (2.2)  are 


Mi  = 


K,  ^  U/ 


(A.8) 


=  (-t>‘ 


(A.9) 


=  (*  ■  I"' 


(A.  10) 
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A  plane  wave  of  angular  frequency  a  is  incident  upon  this  sediment 
configuration  from  above  with  grazing  angle  6.  The  component  of  the  wave 
number  parallel  to  the  fluid/sediment  interface  is 


k 


cos(d) 


(A.11) 


where  pf  and  Kf  are  the  fluid  mass  density  and  fluid  bulk  modulus,  as  given  in 
Table  2.1.  The  simulations  of  the  present  study  concern  the  effectively  one- 
dimensional  case  of  normal  incidence,  with  k  =  0. 

For  the  iit  layer,  the  virtual  mass  coefficient,  c,  and  the  permeability,  re,  of 
Eq.  (2.2)  are  given  by 


c,.  =  1  +  0.5^-^  ,  (A.12) 

and 

<c,  =  M  ,  (A.13) 

where  p/  is  the  pore  size  parameter,  given  by 

ft  =  '  =  1.2,3 . N  .  (A.14) 


3(1  ”  ^N+l) 


(A.  15) 


and  (a)  is  the  mean  of  the  grain  size  distribution. 

The  dynamic  correction  term  F,  of  Eq.  (2.2)  depends  on  the  incident 
frequency.  For  the  iit  layer,  it  can  be  written  as 
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(A.  16) 


Fi 


4 


1  - 


2r(^0 

J 


> 


where 


heFj^d  +  ibeU^d 

ber(li)  +  /bei(^i) 


and 


(A.  17) 


(A.18) 


The  horizontal  position  coordinate  is  denoted  by  x.  The  vertical  position 
coordinate  increases  downward,  and  for  the  iHi  layer  is  denoted  by  z,.  At  the 
upper  boundary  of  the  itti  Biot  layer,  z/=  0.  The  first  vertical  coordinate  zq  within 
the  fluid  medium  is  set  equal  to  zero  at  the  fluid/sediment  interface. 

The  solid  displacement  term  u  and  the  relative  solid/fluid  displacement 
term  w,  from  Eqs.  (2.1)  and  (2.2),  can  be  expressed  in  terms  of  scalar  and 
vector  potentials: 


u  =  .  (A.19) 

w  =  .  (A.20) 

The  scalar  and  vector  potentials  for  the  itb  layer  can  be  written  as 

.  (A.21) 

0^,  =  ,  (A.22) 

~ 3  .  (A-23) 

'ff  =  r/z, j  .  (A.24) 
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Inserting  these  into  Biot's  equations  of  motion,  Eqs.  (2.1)  and  (2.2),  yields  a  set 
of  four  coupled  equations  for  (f>  and  \ia 


where 


o 

II 

1 

+ 

(A.25) 

o 

II 

1 

+ 

(A.26) 

¥si  +  ll  ¥si  =  0  . 

(A.27) 

Vfi  =  YiVsi  . 

(A.28) 

.  F*r]  m 

a.  =  i  -i--'—  , 

(A.29) 

p.  = 

ft  • 

(A.30) 

,,  -  p/®’ 

p,co'  +  a, 

(A.31) 

d,  =  M,//,  -  Cf  , 

(A.32) 

, 

di 

{A.33) 

_  (H.p,-C,p,) 

"  "  ft  “  ' 

(A.34) 

PI  =  ~  , 

(A.35) 

.2  _  (^iPf  ~  ^iPi)  2 

^si  -  ,  1 

d; 

(A.36) 
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= 


(«.p/  -  c,a)  , 

d, 

(A.37) 

2  ,2 

(A.38) 

,2  C,«. 

d. 

(A.39) 

11 

1 

(A.40) 

1 

il 

(A-41) 

The  solutions  to  Eqs.  {A.25)  -  (A.27)  are  plane  waves  of  the  form 


where 


(A.  42) 

+  SuA.,e-“^‘^‘  + 

.  (A.43) 

(A.44) 

+  yl)  -  -  ylif  +  .  (a.45) 


r 

hi 


=  i(4  +  v=)  H-  -  vjf  + 


(A.  46) 


—  ^si  hi 
^si 


(A.47) 


^2i  ”  ^2 

^51 


(A.48) 
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The  terms  Z?/,  (g/,  and  Zs/  are  the  wave  numbers  of  the  Biot  fast,  slow,  and  shear 
waves  within  the  itli  layer,  respectively.  The  "+“  and  coefficients  are  the 
amplitudes  of  the  downward-moving  (transmitted)  waves  and  the  upward- 
moving  (reflected)  waves  within  the  itii  layer,  respectively. 

Let  R  be  the  reflection  amplitude  of  the  wave  within  the  fluid  medium. 
The  unknowns  of  the  problem  are  R,  and  the  As,  Bs,  and  Cs  of 
Eqs.  (A.42)  -  (A.44).  Since  there  are  no  reflections  within  the  homogeneous 
semi-infinite  Biot  medium,  we  have 

~  ^-sN+\  ~  ^-sN+l  ~  ■  (A.  49) 

Thus,  there  are  a  total  of  6ZV  +  4  unknowns,  which  satisfy  the  6AZ  +  4  linear 
algebraic  equations  generated  by  the  boundary  conditions. 

The  boundary  conditions  pertain  to  continuity  of  fluid  pressure,  traction, 
fluid  displacement,  and  solid  displacement.  Within  the  Biot  media,  the  fluid 
(pore)  pressure  is 

ft  =  -  C,(£  -  i  =  i,2 . JV  +  1  ,  {A.50) 


and  within  the  fluid  medium  it  is 

Po  =  p^60'(l  + .  (A.51) 

The  shear  traction  within  the  Biot  media  is  given  by 

(o;L  =  ''  =  1.2 . JV  +  l  ,  (A.52) 

and  within  the  fluid  medium  it  is 

(oi))»  =  0  .  (A.53) 

The  normal  traction  is 
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(A.54) 


(O',).  =  t«,W  -  - ‘'tfe) 

-  C,  (0;  -  e  ^;,)] .  /  =  1, 2, ,..,  W  + 1.  , 


«7.).  =  -A  =  .  (ASS) 

and  the  normal  fluid  displacement  is  given  by 

-  (<Pfi  +  ^ 

i  =  1, 2, N  +  ly 

M,-W,  =  lid -man(ff) .  (A.57) 

Finally,  within  the  Biot  medium,  the  normal  solid  displacement  is  given  by 

(uX  =  (.<!>.  ,  {A.58) 

and  the  tangential  solid  displacement  is 

(uX  =  •  (A-59) 


Neither  of  these  last  two  quantities  have  counterparts  within  the  fluid  medium. 
With  ho  =  0,  all  /V+  1  interfaces  are  subject  to  the  boundary  conditions 


ftUo 

"  ^'-*1  -  A  ’ 

—  n,_i 

i  =  1, 2, ... 

,N  +  1  , 

(A.60) 

/  =  1,2.. 

..,iV  +  l  , 

(A.61) 

“  (^M)2rL  -  h  ’ 

X  "  iz,_j  - 

/  =  !.  2... 

.,iV  +  l  . 

(A.62) 

i 

-  (>^i)=1.0 

=  [(Ui-X  -  (y^i-iX 

lL=c.-’ 

1,2,...,N  +  1  , 

(A.63) 
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and  the  Blot/Biot  interfaces  are  subject  to  the  additional  two  boundary 
conditions  given  by 


L,  r:  0  ^  t  L,..,  =  A,.,  ’ 

/  —  2, 3, . . . ,  1  , 

(A.64) 

i  =  2,  3,  + 1 

(A.65) 

The  transmitted  and  reflected  fast,  slow,  and  shear  wave  amplitudes  within  each 
medium  are  determined  by  substituting  Eqs.  (A. 42),  (A.43),  (A.44),  and  (A.49) 
into  the  above  six  boundary  conditions  and  solving  the  resulting  set  of  linear 
algebraic  equations  for  R,  and  the  As,  Ss,  and  Cs. 
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APPENDIX  B 

REDUCING  THE  COMPUTATIONAL  COMPLEXITY 
OF  THE  PROBLEM 
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PREVIEW 


In  the  process  of  this  investigation,  a  method  was  discovered  that 
dramatically  simplifies  the  computational  complexity  of  the  problem  described  in 
Appendix  A.  Prior  to  this  discovery,  an  A/-layer  problem  required  solving  a 
(6A/+4)x(6A/+4)  banded  matrix  equation.  Problems  with  several  hundred  layers 
would  require  a  few  hours  to  run  on  a  Sun  2000.  With  the  new  simplifications  a 
500-layer  calculation  can  now  be  performed  in  just  over  6  minutes.  This  leap  in 
computational  efficiency  was  made  possible  by  reducing  the  problem  to  a  4x4 
matrix  calculation: 


(MIPSnC)x  =  b  ,  (B.1) 

where  the  4x4  matrix  (M  I  PSnC)  consists  of  a  4x1  matrix  M,  which  depends 
solely  on  the  properties  of  the  fluid  half-space  overlayer,  augmented  with  a  4x3 
matrix  PSnC,  which  is  the  product  of  three  matrices: 

(1 )  a  4x6  prefix  matrix  P  that  is  independent  of  fluid,  sediment,  or 
semi-infinite  Biot  layer  properties, 

(2)  a  6x6  matrix  Sn  that  depends  solely  on  the  A/-layered  sediment 
properties, 

(3)  a  6x3  matrix  C  that  depends  solely  on  the  semi-infinite  Biot 
layer  properties. 

The  4-vector  x  contains  the  solution  for  the  reflection  coefficient  and  the  three 
transmission  coefficients  for  the  Biot  fast,  slow,  and  shear  waves  through  the  N 
sediment  layers.  The  4-vector  b  depends  solely  on  the  properties  of  the  fluid 
overlayer. 

The  sediment  matrix  Sn  contains  all  the  layered  sediment  information, 
and  has  the  following  simpie  form: 
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where  the  6x6  matrix  Lj  depends  solely  on  the  properties  of  the  i^  sediment 
layer.  And  Lj,  in  turn,  can  be  written  as 

L,  =  .  (B.3) 

where  the  6x6  matrix  Bj  depends  solely  on  the  material  properties  of  the  i^^  Biot 
sediment  layer  and  Dj  is  a  diagonal  6x6  matrix  depending  on  both  the  material 

properties  and  the  thickness  of  the  layer.  All  the  layer  thickness  information 
is  contained  in  the  Dj's. 


METHOD 

Given  the  definition  of  the  reflection  and  transmission  coefficients  for  the 
ilh  layer,  as  displayed  in  Eqs.  (A. 42)  -  (A. 44),  we  define  the  vectors  V|  by 


V.  = 


^+si 


i  =  1, 2, 3, ....  N 


(B.4) 


'N+l 


B+sN+1 

yf'+sN+i ) 


(B.5) 


The  values  of  Vj  and  Vj+i  are  related  by  the  boundary  conditions  supplied  in 
Eqs.  (A.60)  -  (A.65).  These  equations  can  be  recast  in  matrix-vector  form  as 


A|V,  =  B,^iV,+,,  i  =  l,2,3,...,N-l  , 


(B.6) 


^  '^N+l  > 


(B.7) 
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where  Ai  and  Bj+i  are  6x6  matrices  describing  the  boundary  conditions  at  the 
lower  surface  of  the  i^^^  layer  and  the  upper  surface  of  the  (i+1)^  layer, 
respectively,  and  C  is  a  6x3  matrix  describing  the  boundary  conditions  at  the 
surface  of  the  Biot  half-space. 

Applying  the  coordinate  system  described  in  the  paragraph  following 
Eq.  (A.18),  the  matrices  for  the  i^i^  layer,  Aj  and  Bj,  with  layer  thickness  hj,  are 
related  by 


where 


A|  —  B,  D, 


D.  = 


fgi'ulU 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


0 

0 

il^h 

0 

0 

0 


0 

0 

0 

g-i‘u  l>i 

0 

0 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


) 


(B.8) 


(B.9) 


and  In,  l2i,  and  Isi  are  the  wave  numbers  of  the  Biot  fast,  slow,  and  shear  waves, 
respectively.  Inserting  Eq.  (B.8)  into  Eq.  (B.6),  and  rearranging  terms,  yields 

V,  =  i  =  l,2,3....,N-l  .  (B.10) 


Applying  this  recursively,  we  have 


V,  =  Dr‘ 


j«»2 


BnV, 


N  ▼N 


(B.11) 


Inserting  Eqs.  (B.7)  and  (B.8)  gives 


Vi  = 


(B.12) 
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Consider  the  matrix 


n  0  0  0  0  0^ 
0  1  0  0  0  0 
0  0  1  0  0  0 
^0  0  0  1  0  0; 


(B.13) 


Note  that  right-multiplying  this  matrix  by  any  6x6  matrix  truncates  the  bottom  two 
rows  of  the  6x6  matrix.  Thus,  if  the  rows  of  Bi,  each  of  which  corresponds  to 
one  of  Eqs.  (A.60)  -  (A.65),  are  ordered  in  the  same  sequence  as  those 
equations,  then  the  4x6  matrix  PBi  describes  the  four  boundary  conditions  at 
the  fluid/Biot  interface  of  the  first  sediment  layer.  Hence,  the  boundary 
conditions  at  the  fluid/Biot  interface  can  be  written  as 

/?M-hPBiVi  =  b  ,  (B.14) 


where  R  is  the  reflection  coefficient,  M  is  a  4x1  matrix  depending  solely  on  the 
fluid  properties,  and  b  is  a  four-vector  which  also  depends  solely  on  the  fluid 
properties.  Inserting  Eq.  (B.12)  into  Eq.  (B.13)  gives 

iJM  +  PS^CvN^,  =  b  ,  (B.15) 


where 


Sn  =  n(B,D;-B;‘) 

J=l 


(B.16) 


Note  that  Eq.  (B.15)  can  be  rewritten  as 

(MIPSnC)x  =  b  ,  (B.17) 

where  the  4x4  matrix  (M  I  PSnC)  consists  of  the  4x1  matrix  M  augmented  with 
the  4x3  matrix  PSnC,  and 
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(  R  \ 


■'^+iW+l 

•®+iN+l 


(B.18) 


Thus,  the  problem  of  Appendix  A  can  be  solved  quite  simply  and  efficiently  by 
forming  the  4x4  matrix  (M  I  PSnC)  and  then  solving  for  x  in  Eq.  (B.19).  When 
compared  to  the  original  (6A/+4)x(6A/+4)  banded  matrix  algorithm,  with 
N  =  500,  this  method  consistently  produced  the  same  results  to  14  out  of  16 
significant  figures,  but  with  a  time  reduction  of  nearly  two  orders  of  magnitude. 
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APPENDIX  C 

CONVERSION  FROM  ONE-  TO  THREE-DIMENSIONAL 
SCATTERING  STRENGTH 


This  page  intentionally  left  blank. 
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The  experimentally  measured  scattering  at  the  water-sediment  interface, 
such  as  the  results  shown  by  Nolle  et  al.,  is  a  three-dimensional  process,  in 
which  energy  is  scattered  within  a  hemispherical  sector  of  2n  steradians,  and 
the  scattered  intensity  decreases  with  distance  due  to  spreading  losses.  The 
simulations  performed  in  this  study  are  one-dimensional,  where  the  scattered 
energy  is  directed  back  to  the  source  of  the  incident  wave  and  there  is  no 
spreading  loss.  In  order  to  make  quantitative  comparisons  between  simulation 
and  experiment,  it  is  necessary  to  derive  a  relationship  between  the  two 
processes. 

The  output  of  the  one-dimensional  simulation  is  essentially  the  reflected 
plane  wave  signal  from  the  sediment  interface.  The  coherent  part  is  identically 
the  reflection  coefficient,  and  the  random  part  may  be  considered  as  the 
scattered  component.  They  are  approximately  equal  to  the  reflection  coefficient 
at  normal  incidence,  that  can  be  experimentally  measured  using  a  real  source 
and  receiver,  after  properly  accounting  for  the  spreading  loss.  The  reflection 
equation  is  as  follows:  Let  Sl  be  the  source  level  of  an  acoustic  source  incident 
normally  on  a  sediment  surface.  Let  Rlc  and  Ru  be  the  coherent  and  random 
parts  of  the  reflected  signal.  Let  both  source  and  receiver  be  at  a  distance  rfrom 
the  sediment  surface.  The  coherent  Rc  and  random  Rr  parts  of  the  reflection 
coefficient  amplitude,  in  decibels,  are  given  by 

Rc  —  Rlc  ~  Si  +20log(2/)+  2cx/'  ,  (C.1) 

Rr  =  Rir  -  Sl  +20  log(2/)  +  2ar  ,  (C.2) 

where  a  is  the  absorption  coefficient. 

When  computed  as  backscattering  strength,  the  corresponding  coherent 
and  random  components  are  defined  as  follows. 

Sc=  Rlc  -  Sl  +40  \og{i)  +  2ar  - 10  log(y4)  ,  (C.3) 

Sr=  Rir  -  Si  +40  log(A)  +  2ar  - 10  log(>A)  ,  (C.4) 

where  A  is  the  insonified  area. 
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Subtracting  Eqs.(C.I)  and  (C.2)  from  (C.3)  and  (C.4),  the  difference  between 
the  reflected  levels  generated  by  the  one-dimensional  simulation,  and  the 
scattered  strengths  measured  in  the  laboratory,  is  given  by 

Rc  -  Sc  =  Rr  -  Sr  =  20  log(2)  +  10  log(>4)  -  20  log{/)  .  (C.5) 

The  insonified  area  is  defined  by  the  beamwidth  of  the  combined  projector- 
receiver  combination.  Let  us  define  the  beamwidth  0  such  that 

/\  =  (0/)2  .  (C.6) 

Then,  the  conversion  from  one-  to  three-dimensional  scattering  strength  is 
expressible  in  terms  of  the  beamwidth  only, 

Rc  -  Sc  =  Rr  -  Sr  =  20  log (2)  +  10  log (6^)  .  (C.7) 

From  the  measurement  of  reflection  loss  (-11  dB)  and  the  peak  of 

backscattering  strength  plot  in  Fig.  23  of  Nolle  et  al.,  the  value  of  Rc  -  Sc  for 

their  experimental  apparatus,  at  1  MHz,  is  found  to  be 

Rc  ■  Sc  =  -11  -22  =  -33dB  .  (C.8) 

Substituting  into  Eq.  (C.7),  the  effective  system  half-beamwidth  at  1  MHz  is 
found  to  be 


6  =  0.011  radians  .  (C.9) 

An  approximate  expression  for  the  combined  -3  dB  beamwidth  of  source  and 
receiver  is  given  as  follows. 

sin(0)V2  =0.443  X/D  ,  (C.10) 

where  X  is  the  acoustic  wavelength  and  D  the  diameter  of  the  aperture,  and  the 
</2  factor  accounts  for  the  cumulative  effect  of  the  projector  and  receiver 

apertures. 
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As  a  check,  substituting  from  Eq.  (C.9)  into  (C.10),  it  is  found  that  the  aperture  of 
both  source  and  receiver  must  have  been  approximately  4  cm,  which  is 
consistent  with  the  stated  transducer  diameter  of  1.5  in.  used  in  the  experiment. 

For  narrow  beams,  the  beamwidth  is  inversely  proportional  to  the 
frequency:  therefore,  the  conversion  from  one-  to  three-dimensional  scattering 
strength  at  500  kHz  is  estimated  to  be 

/=?c  -  Sc  =  -33  dB  +  10  log(4)  =  -27  dB  . 
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